function [cell, gases] = GetNMRCell(temperature, pN2, pXe129, pXe131)
    rb85=Miscellaneous.PhysObject.GetRb85();
    rb87=Miscellaneous.PhysObject.GetRb87();
    n2 = Miscellaneous.PhysObject.GetN2();
    
    xe129 = Miscellaneous.PhysObject.GetXe129();
    xe131 = Miscellaneous.PhysObject.GetXe131();

    T_kelvin=273.15+temperature;
    gases={  Gas.Vapor(rb85, 'temperature', T_kelvin), ...
             Gas.Vapor(rb87, 'temperature', T_kelvin), ...
             Gas.Gas(n2,     'temperature', T_kelvin, 'pressure',  pN2*Torr2Pa), ...
             Gas.Gas(xe129,  'temperature', T_kelvin, 'pressure',  pXe129*Torr2Pa), ...
             Gas.Gas(xe131,  'temperature', T_kelvin, 'pressure',  pXe131*Torr2Pa)};
    cell = CellSystem.NMRCell(gases);
end

